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INTRODUCTION 

The  concept  of  fracture  as  a  result  of  crack  instability  was  intro¬ 
duced  by  Griffith,  who  suggested  two  principles  for  determining  the  on¬ 
set  of  catastrophic  crack  growth  11,2],  One  was  that  fracture  occurs 
when  the  strain  energy  released  by  crack  extension  is  equal  to  the  energy 
absorbed  by  the  resulting  free  surfaces.  The  other  was  that  fracture 
occurs  when  the  tensile  stress  at  the  most  highly  loaded  point  on  a 
crack  surface  reaches  a  critical  value.  The  majority  of  criteria  in 
use  today  for  mixed  mode  fracture  are  related  to  one  or  the  other  of 
these  two  principles.  One  exception  is  a  crude  but  convenient  approxi¬ 
mation  which  assumes  that  only  the  component  of  stress  normal  to  the  crack 
plane  contributes  to  fracture.  Others  include  the  relatively  recent  energy 
density  theory  of  Sih  [3]  and  the  equivalent  specific  energy  theory  of 
Glllemot  [4]. 


A  widely  used  experimental  test  of  fracture  criteria  is  to  measure 
the  strength  and  direction  of  extension  of  a  known  crack  subjected  to 
simple  tension  as  a  function  of  its  orientation.  Unfortunately,  the  pre¬ 
dictions  of  several  criteria  are  within  the  rather  large  egpqriimmral 
scatter,  making  a  choice  difficult. 
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Another  type  of  experiment  that  can  be  used  to  compare  fracture  cri¬ 
teria  employs  unknown,  naturally  occurring  cracks,  and  Is  statistical  in 
nature.  A  major  objective  of  statistical  theories  of  fracture  is  to  use 
the  experimentally  established  fracture  behavior  in  simple  tension  to  pre¬ 
dict  fracture  behavior  in  some  other  stress  state.  Since  the  result  de¬ 
pends  on  the  fracture  criterion  assumed  [5],  comparison  of  the  prediction 
with  experiment  can  be  used  to  eyaluate  the  relative  merits  of  fracture 
criteria.  A  convenient  second  stress  state  for  this  purpose  is  equibi- 
axlal  tension.  The  purpose  of  the  present  paper  is  to  carry  out  such  a 
comparison. 

THEORY 

Assume  that  failure  can  occur  due  to  any  one  of  many  independent  and 
mutually  exclusive  mechanisms  or  causes,  each  having  the  infinitesimal 
probability  of  failure  (APf)^.  Under  these  circumstances,  the  probability 
of  survival  can  be  shown  to  be  [5] 

Ps  «  exp[-Z(APf)i]  (1) 

In  the  case  of  a  body  containing  Griffith  cracks  (open  or  closed  ellip¬ 
tical  cylinders)  randomly  oriented  about  axes  normal  to  the  stress  plane 
and  subjected  to  a  uniform  state  of  stress  E 
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Here  V  *  volume 

o  *  critical  stress,  defined  as  the  remote  tensile  stress  that  will 
cause  fracture  when  applied  normal  to  the  crack  plane. 

N(a  )  =  number  of  cracks  per  unit  volume  with  critical  stress  <  0c« 
w(E,a^)  =*  angle  within  which  normal  to  crack  plane  must  lie  to  cause  fracture. 

The  right  hand  side  of  (2)  is  simply  the  probability  that  a  crack  in  the 
critical  stress  range  Aoc  is  present  times  the  fraction  of  such  cracks  with 
an  orientation  that  will  result  in  fracture.  Inserting  (2)  into  (1)  and 
going  to  the  limit  of  small  Ao 


Pg(E)  -  exp 


We  note  that  the  integral  in  (3)  is  finite  because  <o  *  0  when  0C  >  o i, 
the  maximum  principal  stress  (in  some  cases,  to  be  discussed  later,  0 
may  exceed  0^  by  several  percent  for  a>  to  be  zero) . 

In  the  case  of  equibiaxial  tension 

w  *  ir  (Oc  <  a  )  (4a) 

=  0  (©c  >  0  )  (4b) 

As  a  result,  (3)  reduces  to  the  simple  form 

Pg(o,0)  =  exp[-VN(o) ]  (5) 

To  facilitate  the  comparison  between  uniaxial  and  equibiaxial  fracture 
behavior,  we  assume  with  Welbull  [6]  that  N  takes  the  form 


N(0c) 


kO 


m 


Then 


where 


Pg(0,0)  =  exp  jj-Vkom  J 

Pq(o,o)  «  exp  Jj-Vkm  Jq  0  ®_1  dOcj  *  exp^-Vk’o®  ] 
m-l 

J  o  ir  x 


k  -  km 
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(7) 

(8) 

(9) 

(10) 


The  desired  relation  between  uniaxial  and  equibiaxial  fracture  statistics 
can  thus  be  expressed  in  the  form 


£nPg(0,0) 

JtnPg(0,o) 


k  f  f00  0)  m-l  1-1 

k’"[mJo  ¥x  dxJ 


(11) 


As  a  simple  illustration  of  the  use  of  (11),  consider  the  case  of 
shear-insensitive  cracks.  With  such  cracks,  fracture  depends  only  on 
the  normal  component  of  stress  0  which  for  simple  tension  is  given  by 


0  ■  0sin  $ 

A 


(12) 
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Figure  2.  Uniaxial/equibiaxial  relation  for  five  fracture  criteria. 


Other  fracture  criteria  leading  to  other  functional  relations  between 
w  and  x  and  therefore  other  functions  Fj.(m)  have  been  discussed  in  [5]  and 
(81 •  For  the  sake  of  brevity,  we  content  ourselves  here  with  the  identi¬ 
fication  of  crack  types  and  fracture  criteria  leading  to  curves  F~  to  Fc, 
as  shown  in  Table  I. 


CURVES 

CRACK  TYPE 

FRACTURE  CRITERION 

F2(m) 

Open;  Griffith 

Tensile  stress  in  crack  sur¬ 
face. 

Fj  (m) 

Open;  Penny-shaped 

F4(m) 

Open  or  closed; 

Griffith 

Energy  release  rate  for  co- 
planar  extension. 

F5(m) 

Open  or  closed; 
Penny-shaped 

Table  I.  Identification  of  curves  in  Fig.  2. 


The  curves  for  F,  and  F.  depend  on  Poisson's  ratio  v,  which  was  assumed 
to  be  0.25.  i  3 
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In  addition  to  to  F5,  Fig.  2  shows  experimental  results  for  three 
different  materials.  The  circle  represents  test  results  taken  on  tubes  of 
pyrex  glass  [9],  The  triangle  represents  test  results  on  ATJS  graphite 
bars  and  disks  spun  to  failure  [10],  In  the  case  of  the  disks,  only  the 
failures  originating  very  near  the  center  were  included  in  the  analysis. 
The  square  represents  results  on  beams  and  disks  of  alumina  subjected  to 
bending  stress  [11],  The  test  results  are  summarized  in  Table  I. 


Material 

m 

k/k' 

Pyrex 

3.35 

2.0 

ATJS  Graphite 

12.00 

3.25 

Alumina 

14.3 

3.35 

Table  II.  Experimental  k/k*  relation. 


It  is  evident  that  agreement  with  the  experiment  leaves  something  to 
be  desired.  We  next  turn  to  predictions  based  on  energy  density  theory 
and  of  fracture  criteria  strain  energy  release  rate  for  non-coplanar  crack 
extension. 


Sih  [12]  has  shown  that  for  an  inclined  Griffith  crack  in  simple  tension 
the  critical  energy  density  function  is  given  by 


Sc  “  all\2  +  2a12klk2  +  a22k22 

(16) 

where 

an  ”  T27T  [3-4V-COS0  )  (1+cosG  )] 
xoy  0  0 

(17a) 

a12  =  2sin0o[cos0o-(l-2v)] 

(17b) 

a,,  *  TZTT  I4(l-v)(l-cos0  )+(l+cos0  )(3cos0  -1)] 
xz  lop  0  00 

(17c) 

and  where  y  is  the  shear  modulus.  For  Griffith  cracks 

kx  -  aa1/2sin2e 

(18a) 

1/2 

k£  ”  oa  sing  cos8 

(18b) 

while  for  penny-shaped  cracks 

,  2  „  1/2  .  2~ 
k  *  —  0a  sin  p 

1  TT 

(19a) 

.  4  _  1/2  sin8  cos8 

k,  -  —  0a  — r— - 

L  TT  Z— V 

(19b) 

6q  is  chosen  by  minimizing  S  with  respect  to  9. 
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A  plot  of  the  resulting  x(0)  is  shown  in  Fig.  3.  Also  shown  is  a  plot  of 
the  corresponding  relation  for  penny-shaped  cracks. 


A  major  motivation  underlying  Sih's  early  efforts  in  developing  energy 
density  theory  was  the  fact  that  up  to  that  time,  mathematical  difficulties 
had  precluded  accurate  calculation  of  strain  energy  release  rates  for  non- 
coplanar  crack  extension.  Recently  two  solutions  have  been  offered  for  the 
case  of  closed  Griffith  cracks,  one  by  Palaniswamy  and  Knauss  [14],  and  one 
by  Wu  [IS].  The  two  solutions  are  in  excellent  agreement  with  each  other, 
and  lead  to  virtually  the  same  strength  for  an  angled  crack  as  that  obtained 
from  the  theory  proposed  by  Erdogen  and  Slh  [16]  using  Griffith's  stress  cri¬ 
terion.  This  strength  can  be  readily  calculated  from  the  following  relation  [17] 


a  6  /  0 

~  "  cos  2~  sin0lsing  cos2  -t~ 


-  4  cos0  sin© 

2  o  i 


(19) 


where  0q  is  found  from  the  equation 

tan$  sin9  +  3cos9  =1.  (20) 

o  o 

The  resulting  function  x(B)  =  cfc/a(0)  is  shown  as  the  curve  in  Fig.  3, 
labeled  "E.  and  S."  Note  that  this  curve  reaches  a  maximum  value  of  about 
1.03  at  Bc  ~  70°,  whereas  all  the  others  considered  here  reach  a  maximum 
value  of  x  ■  1  at  8  «  90°.  This  peculiarity  has  important  consequences  for 
the  corresponding  value  of  k/k*. 

When  the  relations  of  Fig.  3  are  used  to  evaluate  to  and  the  results 
substituted  in  (11),  the  curves  for  k/k'  shown  in  Fig.  4  are  obtained.  The 
energy  density  criterion  leads  to. fairly  good  agreement  with  the  experiment 
in  the  case  of  Griffith  cracks,  and  quite  good  agreement  for  penny-shaped 
cracks.  The  agreement  of  statistical  theory  based  on  the  use  of  available 
theories  for  non-coplanar  energy  release  rate  is  so  poor  as  to  be  (in  the 
opinion  of  the  present  author)  unacceptable  because  it  leads  to  an  incorrect 
variation  of  k/k'  with  m. 
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DISCUSSION 

The  preceding  analysis  considers  open  or  closed,  Griffith  or  penny¬ 
shaped  cracks,  fracturing  according  to  any  one  of  nine  different  fracture 
criteria  (Fg  does  double  duty,  representing  both  the  Og  and  the  non-coplanar 
energy  release  rate  criteria).  Best  agreement  with  the  experiment  was  ob¬ 
tained  using  energy  density  theory  as  applied  to  closed  penny-shaped  cracks. 
The  poor  agreement  of  the  relation  Fg  with  the  experiment  is  quite  surprising 
since  it  appears  to  offer  in  principle  an  accurate  mathematical  solu.ion  of 
a  physically  valid  principle  (energy  balance  for  non-coplanar  crack  exten¬ 
sion).  The  low  value  for  the  ratio  k/k*  can  be  traced  to  the  influences 
on  the  integral  in  (9)  of  the  portion  of  the  x(£)  curve  that  is  greater 
than  unity.  This  is  not  only  physically  implausible,  but  leads  to  the 
incredible  prediction  that  as  m  -*■  ",  k/k'  -*■  0. 

It  was  to  be  expected  that  of  the  two  energy  density  solutions,  the 
one  based  on  use  of  penny-shaped  cracks  is  in  best  agreement  with  the 
experiment  since  such  cracks  are  more  closely  related  to  naturally  occurring 
cracks  than  Griffith  cracks.  To  validate  the  conclusion  that  it  agrees  with 
uniaxial/biaxial  strength  data  better  than  any  other  well-known  theory, 
however,  it  is  necessary  to  justify  the  use  of  two-dimensional  statistical 
theory  (eq.  (11))  rather  than  the  three-dimensional  theory  that  might  appear 
more  appropriate  [5],  This  justification  lies  in  the  surprising  fact  that 
the  two  theories  lead  to  the  same  k/k'  ratio.  This  was  proved  in  [16]  for 
the  relations  Fi  and  F^  (see  Fig.  2),  and  has  subsequently  been  shown  by 
the  present  author  to  be  a  generally  valid  result. 

It  should  be  pointed  out  that  the  present  paper  is  not  the  first  to 
discuss  statistical  fracture  theory  using  energy  density  theory  as  the 
fracture  criterion.  In  a  couple  of  papers  published  in  1977  [19,  20] 
Jayatilaka  and  Trustum  considered  uniaxial  and  biaxial  stress  states, 
respectively.  However,  their  results  are  flawed  by  the  use  of  a  crack 
size  distribution,  and  therefore  a  functional  relation  for  N(ac),  that 
was  determined  experimentally  for  surface  cracks  in  glass.  There  is  no 
reason  to  expect  that  the  same  law  should  hold  for  interior  cracks  in 
polycrystals.  In  fact,  it  is  very  likely  it  would  not  hold  even  for 
another  sample  of  the  same  glass,  since  surface  flaws  in  glass  are  gen¬ 
erally  considered  to  be  due  to  handling  damage.  In  addition,  no  com¬ 
parison  of  their  results  were  made  with  experimental  data,  or  even  with 
predictions  of  statistical  fracture  behavior  based  upon  other  fracture 
criteria. 

Finally,  there  is  the  question  of  the  reliability  of  the  data.  The 
fact  that  the  three  data  points  (a)  each  represent  many  tests,  (b)  in¬ 
volve  different  types  of  experiments,  (c)  were  carried  out  by  different 
experimenters  on  different  materials,  and  (d)  collectively  establish  a 
trend  closely  related  to  what  would  be  expected  on  theoretical  grounds, 
tends  to  establish  confidence  in  their  validity.  It  would  nevertheless 
be  highly  desirable  to  obtain  additional  experimental  data  to  give  greater 
assurance  to  the  conclusions  reached  in  this  investigation. 
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C0NCLUSI0NS 

1.  The  relation  between  the  fracture  statistical  behavior  of  specimens 
in  simple  tension  and  that  in  equibiaxial  tension  depends  on  both 
the  assumed  crack  shape  and  the  fracture  criterion  employed. 

2.  The  present  study  investigates  the  uniaxial/ equibiaxial  relation  for 
Griffith  and  penny-shaped  cracks  for  a  number  of  fracture  criteria. 
Best  (in  fact,  very  good)  agreement  with  the  experiment  was  obtained 
assuming  energy  density  (or  absorbed  specific  energy)  theory  and 
penny-shaped  cracks. 

3.  Contrary  to  what  might  be  expected,  the  maximum  energy  release  rate 
criterion  for  non-coplanar  extension  of  Griffith  cracks  leads  to  a 
uniaxial/equibiaxial  fracture  relation  that  is  in  poor  agreement 
with  experiment. 
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